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© System for the rapid convergence of an adaptive filter In the generation of a time vacant signal for 
cancellation of a primary signal. 

© System for generating a secondary signal (y s (k|) for suppression of a primary signal (y p (k)), Comprising 

- means (8) for generating a reference signal 

- a control unit (1) having a digital filter for providing a cancellation control signal (s(k)); 

- cancellation-generating means (2) for generating a secondary signal (y s (k)); 

- sensor means (4) for providing an error signal (ffk)^ 

- update means (5) for providing an update signal (up(AQ); 

- an adaptation unit (7) for supplying a filtered-x signal [r{k)) to the update means (5); 

- means for updating filter coefficients {e c .Jk)) of the filter in the control unit (i) with the aic^ of the so-called 
"projection algorithm" and 

- the various features being such that an input correlation matrix defined according to equation 17 has an 
eigenvalue distribution which has a value which is at least substantially identical to 1 aftei as few steps as 
possible. 
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The invention relates to a system for the generation of an estimated time variant signal for suppression 
of a primary signal, comprising: 

- means for generating at least one reference signal; 

- a control unit at least provided with a digital filter, a first input for receiving the reference signal, a 
5 second input for receiving an update signal for updating coefficients of the digital filter and an output 

for providing a cancellation control signal; 

- cancellation-generating means which are connected to the output of the control unit for the generation 
of a cancellation signal, which is intended, after propagation along a secondary transfer path having a 
secondary path transfer function, to be added as the estimated time variant signal at an addition point 

w to the primary signal in order to provide a residual signal, 

- sensor means for measuring the residual signal at the addition point and for providing an error signal; 

- update means provided with an input which is connected to the sensor means and an output for 
providing the update signal. 

A system of this type is disclosed in US Patent 4 878 188. 

is The known system can be used, inter alia, to cancel out irritant noise (the primary signal) from a rotary 
machine and comprises means for generating a reference signal, which means, for example, comprise a 
tachometer of a rotary machine or an optical sensor for measuring the speed of revolution of the rotary 
machine. In addition, the known system comprises a signal processor for generating a cancellation control 
signal, which is supplied to cancellation signal-generating means for generating a cancellation signal, which 

20 as far as possible has to cancel out the primary signal within a specific area. The residual signal in that area 
is measured using suitable sensor means, the output signal from which is fed to update means. Said latter 
means produce an update signal for the signal-processing means in order to update the filter coefficients of 
an adaptive filter present therein. The update signal is derived form the residual signal and can, for 
example, be determined with the aid of the least mean squares (LMS) method and is aimed at causing the 

25 residual signal to become as small as possible. 

In the system according to the said US Patent, in the case of a multi-channel embodiment, the 
reference signal is supplied via n base reference signals on n parallel branches to the signal-processing 
means. Each base reference signal is fed to a cosine generator and a sine generator. The cosine generator 
produces a base cosine signal with a specific base frequency and a number of higher harmonics thereof. 

30 The sine generator produces a base sine signal with a specific base frequency and a number of higher 
harmonics thereof. All cosine and sine signals are multiplied by suitable filter coefficients of the adaptive 
filter. Phase means are used to ensure that the adaptation of the filter coefficients takes place in such a way 
that the filter characteristics remain within 90 • of the primary signal. A summation unit finally adds together 
all cosine and sine signals modified in this way in order to generate the cancellation control signal. 

35 The convergence speed of the algorithms used in the prior art is found to be moderate for various 
types of primary signals, for example for rapidly changing periodic signals which, for example, arise in the 
exhaust of a car, and for pulse-like signals. In this context the convergence speed is defined as the speed 
with which a control system used proceeds to the state of as complete as possible noise suppression. The 
convergence speed of the known algorithms is found not to be optimum in multi-channel systems as weil. 

w One objective of the invention is to provide a system and a method with which the speed of 
convergence of the algorithm used is increased compared with the prior art. 

To this end, a system of the type mentioned in the preamble is characterised in that 

- the system also comprises an adaptation unit which has an input for receiving the at least one 
reference signal and an output for providing a filtered-x signal, which output is connected to a further 

45 input of the update means; 

- the adaptation unit has a transfer function which is the same as the transfer function of the secondary 
path; 

- means are provided for updating the fitter coefficients with the aid of the so-called "projection 
algorithm" and 

so - the means for generating the reference signal generate the said reference signal in such way that an 
input correlation matrix defined as follows 

y flk)T*W 

55 hi T*{k) 7(Jt) 

has an eigenvalue distribution which has a value which is substantially equal to 1 after as few steps as 
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possible. 

As a result of the use of a combination of a projection algorithm and the selection 
signal in the manner described above, in practice a more rapid convergence is achieved 
according to the prior art. 

Preferably, the system according to the invention is a multi-channel system, that is to 
which the cancellation-generating means comprise various cancellation-generating units 
means comprise various sensor units, the control unit being provided with means for 
of the digital fitter in accordance with the following recursive relationship: 
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where: 



15 is an L-dimensional filter coefficient vector when L cancellation signal generators are present; 

7 = step size parameter; 

R„(k) is an LxM matrix which contains input signals r v ,^ m (k) as elements when M sensor^ are present; 
t(k) = error signal vector with M elements originating from the M sensors 

In an embodiment for a single channel system, the system according to the invention is characterised 
so in that 

the means for generating the reference signal are equipped to generate the following j-eference signal: 
25 where: 

- the update unit is provided with means for calculating the update signal as follows: 



UPy{k) 



i v {k)e(k) 



U' + rUk) r v {k)) 



where: 

7 = step size 



a' = small constant to prevent division by 0; 
- the control unit is provided with means for calculating the filter coefficients as follows: 



QcM = e^k-^-up^k) 



The system according to the invention can also comprise estimating means for estimating 
function of the secondary path, which means have an input connected to the sensor means 
connected to the adaptation unit for transmitting an output signal for updating of the 
in the adaptation unit. 

The control unit, the update means and the adaptation unit, and the estimating 
optionally used, can be implemented together in the software of a computer. 

The cancellation-generating means can be selected form the following means: 

at least one loudspeaker, at least one vibration generator, at least one adjustable 
source or a combination thereof, and the sensor means can be selected from the following 

at least one noise pressure sensor, at least one vibration recorder or a combination thereof. 

The invention will be explained in more detail below with reference to a few 
which are intended to illustrate the principle of the invention and not to imply any restrictior 
Figure 1 shows a block diagram of a known anti-noise system; 
Figure 2 shows an example of a block diagram of an anti-noise system according to the 

The following is known from mathematics- Assume there is a primary function P(n): 
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Pin) = gai<") -filxj (1) 

5 

where: 

ajfn) = multiplication coefficient 
Fi(x„) = a known function having a variable x„ 
I = 1, 2, 3, I where I has a predetermined value. 
jo A number of conditions must be met if it is desired to calculate the various aj(n) from a number of 
measurements of x n . The first condition is that: I S n and the second that the various ai(n) must be 
stationary. The various a,(n) can then be calculated using the least squares method. The number of steps 
needed for this is equal to the number of coefficients a^n). 

However, if P is contaminated with noise, the least squares method can then not be used. Provided the 
is other two abovementioned conditions are met, the least squares error method can then be used. 
Furthermore, large matrices can sometimes occur in the calculations, which matrices impose a heavy 
burden on the computing time of a computer used, and the primary function P(n) is far from being always 
available, as in anti-noise systems. If this is the case, the recursive least squares error method could be 
used. If the coefficients a ( (n) are themselves not time-invariant, the least squares method can not be used at 
20 all. 

Should P be contaminated with noise, should a particularly large matrix occur in the calculation and 
should the coefficients a,(n) themselves not be time-invariant, the so-called "projection algorithm" can then 
be used successfully, as is known per se from mathematics. By using a projection algorithm together with 
orthogonal functions, the I coefficients ai(n) can even be estimated in n steps. A direct application of a 

25 projection algorithm with orthogonal functions in a system for the production of, for example, anti-noise is 
not possible without the necessary adjustments because the path traversed by the anti-noise (the 
"secondary path") introduces phase shifts and time delays. Furthermore, in general, orthogonal functions 
are not supplied to the system. The necessary adjustments to a projection algorithm with orthogonal 
functions for, for example, an anti-noise system will be explained below. It will be seen from this explanation 

so that, In practice, the definition of orthogonal functions frequently takes up too much computing time and an 
approximation of the ideal state - that is to say the use of a projection algorithm with orthogonal functions - 
is to be preferred. 

Figure 1 shows a known model of an anti-noise system. A reference signal <f>(k) is generated by a 
reference signal generator 8 and is fed to a control unit 1, which calculates a cancellation control signal s(k). 
35 The cancellation control signal is fed via suitable connecting means to a cancellation signal generator 2, 
which generates a cancellation signal which, after passing over a secondary path with transfer function H 
arrives as process output signal y s (k) in a specific area 3, where it has to cancel out a primary signal y p (k) 
by means of addition. The signal resulting after addition - the residual signal y,(k) - is detected with the aid 
of a sensor 4, which feeds an output signal «<k) to update means 5. The update means 5 also receive the 
40 reference signal <^(k) and produce an update signal up(k), which is fed to the control unit 1, where it is used 
to update the filter coefficients of a filter used. 

In a multi-channel system there are various cancellation signal generators 2 and various sensors 4. 
The primary signal y p (k) can be modelled as a derivative of the reference signal $(k), said latter signal 
traversing a primary path with primary transfer function H p . If it is assumed that the reference signal ^(k) 
45 can be represented with sufficient accuracy by N frequency components, then <£(k) can be represented as a 
single reference signal in the complex plane as follows: 

q> <ic> = f;exp(j2Jtf v JrAC) (2) 

50 v . 3 



where: 

4>(k) = reference signal at time kAt 
55 f u = a frequency component 

The primary signal y p (k) will have traversed the primary path before it has arrived at point 3. If it is now 
possible to use a finite impulse response of length I for a primary process associated with the primary path 
and for the control unit 1, then both the primary signal y p (k) and the cancellation control signal s(k) can be 

4 



BNSDOCID: <EP 



0654901A1,I_> 



EP 0 654 901 AT 



10 



represented as linear combinations of successive signals *(k). *(k-1) ... cMk-l + 1). The input elements can 
be collected in a time vector <Mk): 

Tm = [<M*W>(/(-1)...*(K-/+1)] 7 (3) 

The filter based on a finite impulse response which is used in the control unit 1 is usually termed a 
transversal filter. 

The reference signals can also be defined differently. Specifically, independent wave form generators 
can be used in order to generate various reference signals. The input signals for the prim ary process with 
transfer function H p and the control unit are then multidimensional. If it is assumed that tre primary signal 
comprises I frequency components v = 1....I, then the primary signal y p (k) and the cancellation control 
signal s(k) comprise a linear combination with complex coefficients of the following time functions in the 
complex plane: 

rs HK) = exp(j2vfMt) (4) 

which are the functions on which a Fourier transform can be based. The said complex coefficients are then 
Fourier coefficients. The input elements can thus also be modelled as a "Fourier vector' 

20 UK) = [*,(A)* z (*)...4>i(*)J r (5) 

A control unit which is based on equation (5) can be termed a Fourier control unit. 
If it is to be possible to apply the abovementioned projection algorithm with orthogorjal 
successive reference vectors according to equation (3) or (5) (that is to say for, for 
k=j + 1) will always have to be orthogonal. For an algorithm which is based on reference 
defined according to equation (3) it is virtually never possible to make successive 
orthogonal. In principle, this is possible for two successive Fourier vectors 4>i<A) as define<jl 
However, in practice this takes a very great deal of computing time and would not serve 
convergence. The latter is particularly true if the orthogonalisation Is carried out in the 
least mean squares algorithm (for details on this see: B. Widrow and S.D. Stearns, 
Processing", Englewood Cliffs, Prentice Hall, 1985; S.J. Elliott, I.M. Stothers and P.A. 
error IMS algorithm and its application to the active control of sound and vibration", IEEE 
Speech, Signal Processing, Vol. ASSP 35, pp. 1423-1434, Oct. 1987; and L.J. Erikssor 
R.A. Greiner, "The selection and application of an MR adaptive filter for use in active 
IEEE Trans. Acoust, Speech, Signal Processing, Vol. ASSP 35, pp. 433-437, April 1987) 

The problem of too long a computing time can be avoided by selecting or modifying 
vectors in such a way that the eigenvalue distribution of an input correlation matrix for the 
has become substantially equal to 1 after as few steps as possible, which steps are needed 
filter coefficients of the filter present in the control unit even though this will take place in 
is the case with an orthogonal system of reference vectors. This will be explained in more 

For a series of reference vectors ^(/r) (k = 1, 2, ...), as defined in equation (3 
distribution of this type is usually found to be far from optimum. This is frequently still true 
of steps of the algorithm used. However, it is found that if equation (5) and the series of 
Uk) ( k = 1.2. •••) as defined in equation (4) are used, the eigenvalue distribution of 
matrix is already equal to 1 after 1 signal period. 

The preferred algorithm will now be derived for a single channel anti-noise system 
where both the primary process and the secondary process are linear. 
The cancellation control signal s(k) (k= 1,2 n) can be written as: 



25 



30 



35 



40 



so s(k) = ^ H (k)9^) (6) 

wherein " H " indicates the transposition of the conjugate of the complex vector, ^(k) is 
comprises the input values for the control unit 1 and *c(A-1) is the vector which comprises 
filter coefficients of the control unit 1. 
55 If it is assumed that the time-invariant secondary process can be modelled as a finite i 
(FIR) with coefficients hj, then the output signal y s (k) of the secondary process can 
convolution of the cancellation control signal s(k) and the finite impulse response fy: 
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y e {k) = "PhjSik-j) {7) 



For an ideal sensor 4: the measured error t(k) = y p (k) + y s (k) so that tor a time delay control unit the 
error «(k) can be written as: 



10 



ilk) = y„<M + Yh J V 9 '(k-j-i)6 <! , i lk-j-l) 



(8) 



In this equation denotes the complex conjugate of the vector concerned. For a Fourier control unit based 
15 on equation (5) the error e(k) can be written as: 



"a v-i 



20 



In order to simplify the estimation process, use is made of the abovementioned filtered-x configuration, 
which is known per se. In this context it is assumed that if the step size parameter 7 is chosen to be small, 
the change in the estimated filter coefficients is small over the length of the finite impulse response (over J 
25 samples). In other words, the following equation then applies: 

- M*-1 ) for ;= 0 J-1 (10) 

Based on equation (1 0), the error «(k) can be written as: 

30 - 



35 

when a time delay control unit is used and as: 



(11) 



40 



J-1 



elk) * y p {k) *Ee c , v (M)T;v^-il 

V.1 J^O 



(12) 



when a Fourier control unit is used. Equations (11) and (12) can also be written in the following vector 
45 notation: 

«<ft>«ypM+^*)M*-i) (13) 

^ in which the components of the vector r"(K) are equal to the value defined in equation (14) in the case of a 
so time delay control unit and are equal to the value defined in equation (1 5) in the case of a Fourier control 
unit: 
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I,{k) = £ A/cp (k-j-i) 
r v {k) = ^h'vjk-j) 



(14) 
(15) 



j -a 



The vector Ftyc) is termed the fiitered-x signal. 

If the secondary path is known, updating of the filter coefficients is reduced to a 
configuration if the step size parameters are small. If the projection algorithm, which is 
applied to the abovementioned filtered-x configuration, the following algorithm can be 
of the filter coefficients: 



general linear control 
know n per se, is now 
derived for updating 



W c lk) = ITjk-i) 



r(k)e(k) 
a*z"(k)l{k) 



(16) 



where: 

fl e (K) = filter coefficient vector at "time" k 
7 = step size parameter 

o = a small positive constant to prevent division by 0 

Updating of the filter coefficients with the aid of algorithm (16) can be implemented 
diagram shown in Figure 2. In Figure 2, elements which are the same as those in Figure 
by the same reference_numerals. Instead of a reference signal 4>(k) reference is ma^e 
reference signal vector $(K), as defined in either_ equation (3) or equation (5) 
single update signal up(k) an update vector up{H) is shown in which all update si 
coefficients of the control unit 1 are incorporated. An adaptation unit 7 has been 
reference signal 4>*(/f) and the input of the update unit 5. The adaptation unit 7 adapts the 
<t>"(k) supplied to the update unit 5 in accordance with the processing of equation (14) or " 
whether the reference signal is defined according^ to equation (3) or equation (5) respect 
supplied to the update unit is thejiltered-x signal f(k). On the basis of said filtered-x signa 
unit calculates an update signal up(k) in accordance with the following equation: 



up(k) 



rlk)e(k) 
a +F 5 (Jc)7(ic> 



(16') 



The normalised input correlation matrix in this filtered-x configuration is equal to (k = 1 



z(k)z"i.k) 



(17) 



* r"ik)T{k) 
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It is now found that the eigenvalue distribution of this autocorrelation matrix is still 
the input vectors. Optimisation is generally difficult for a time delay control unit. For a 
however, the matrix according to equation (17) can easily be optimised if the matrix were 
matrix, that is to say if all elements of the matrix not lying on the diagonal were to be 
the diagonal elements. This condition is met if two arbitrary components of the input 
independent. It is probable that for a Fourier control unit the elements of the matrix (1 7) 
diagonal will have become negligibly small with respect to the diagonal elements after a 
calculation, so that in this case the matrix (17) is a good approximation to a diagonal matrb 
distribution of the matrix (17) will then be substantially equal to 1 if the componentse c , v (fr) 
in the update algorithm according to equation (16) are scaled by the following factor: 
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Use of equation (18) results in the following recursive relationship for the calculation of the filter coefficients 
0c.M- 

Q (Jc)=6 Ik 1) t r * lk)eik) (19) 

rs wherein: a' is a small constant to prevent division by 0. 

In the case of implementation of the algorithm according to equation (19), the real component must be 
separated from the imaginary component. If only one filter coefficient 0 c (/c) were to be calculated, algorithm 
(19) would be reduced to: 



(RelQ e (k))) (Reie e (k-1)}) {¥m\r\k)))* {k) (2 0) 

\lm[e c (k)]j ' \lmie c (k-l))j- y a + Relr(k)] z + Imlr[k)] 2 

25 

In the case of multiple input signals algorithm (20) applies for each e c ,»(#(). 

For multi-channel systems, that is to say systems in which diverse cancellation signal generators 2 and 
diverse sensors 4 are used, a similar algorithm can be derived for updating the filter coefficients in the 
control unit i . For a multi-channel system, the following recursive relationship can be derived for the filter 
30 coefficients fl ClV (/f): 

B7;<*) =B^U-l) -ylR Y lk)R?(k))- l R v (k)e(k) < 21) 

35 

where: 

sXv(/r) 'S an L-dimensional filter coefficient vector when L cancellation signal generators 2 are present; 
7 = step size parameter; 

Rv(k) is an LxM matrix which contains input signals r Vitim (k) as elements when M sensors 4 are present; 
40 c(k) = error signal vector with M elements originating from the M sensors 4. 

It can be demonstrated that the convergence speed is also optimum for the multi-channel algorithm 
according to equation (21) and is not dependent on the position of the various cancellation signal generators 
2 and the various sensors 4. 

In the above it has been assumed that the transfer function H of the secondary path is known and 
45 stationary. In reality, the secondary path is far from always being stationary, which means that the 
components hj of the transfer function H, as used in the adaptation unit 7, are not stationary. Means 6 for 
estimation of the transfer function at the time of updating of the filter coefficients can also be incorporated in 
the diagram. Said estimating means 6 can be installed as shown in Figure 2, that is to say connected by an 
input to the sensor 4 and by an output to the adaptation unit 7 for updating of the components hj of the 
so transfer function H. Means of this type are known to those skilled in the art and do not require more detailed 
explanation here (see, for example, US Patent 4 878 1 88). 

It will be clear to a person skilled in the art that the control unit 1. the update unit 5, the estimating unit 
6 (if used) and the adaptation unit 7 can be implemented together in the software of a computer. 

55 Claims 

1, System for the generation of an estimated time variant signal (y s (k)) for suppression of a primary signal 
(y P (k)), comprising: 
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means (8) for generating at least one reference signal (4>(k)); 
a control unit (1) at least provided with a digital filter, a first input for receiving 
signal U>{k)), a second input for receiving an update signal (up{k)) for updating 
(k)) of the digital filter and an output for providing a cancellation control signal 
cancellation-generating means (2) which are connected to the output of the control 
generation of a cancellation signal, which is intended, after propagation 
transfer path having a secondary path transfer function (H), to be added as th^ 
variant signal (y s (k)) at an addition point (3) to the primary signal (y p (k)) in 
residual signal (y,(t)), 

sensor means (4) for measuring the residual signal (y,(t)) at the addition 
providing an error signal («<k)); 

update means (5) provided with an input which is connected to the sensor 
output for providing the update signal (up(k)), 
characterised in that 

the system also comprises an adaptation unit (7) which has an input for receiving 
reference signal (<>(*)) and an output for providing a filtered-x signal (r(/c)), 
connected to a further input of the update means (5); 
the adaptation unit (7) has a transfer function (H) which is the same as the 
the secondary path; 

means are provided for updating the filter coefficients (8 c , v (k)) with .the aid 
"projection algorithm" and 

the means (8) for generating the reference signal (<t>{k)) generate the said 
way that an input correlation matrix defined as follows 



p zlk)7*(k) 
3 ~P(k) x(k) 



means 



the at least one 
which output is 

transfer function (H) of 

of the so-called 

I reference signal in such 



has an eigenvalue distribution which has a value which is substantially equal 
steps as possible. 



System according to Claim 1, characterised in that the cancellation-generating 
various cancellation-generators and the sensor means (4) comprise various sensors 
being provided with means for updating coefficients of the digital filter in accordance 
recursive relationship: 



generators (2) are 
M sensors (4) are 



where: 

e c ,y{k) is an l-dimensional filter coefficient vector when L cancellation signal 
present; 

45 7 = step size parameter; 

R v (k) is an LxM matrix which contains input signals r Vj1>m (k) as elements when 
present; 

*{k) = error signal vector with M elements originating from the M sensors (4). 

so 3. System according to Claim 1 , characterised in that 

- the means for generating the reference signal (<*>(*)) are equipped to generate the following 
reference signal (4(k)): 



where: 

<MK) = exp(/2w/,Mf) 



the reference 
^efficients (9 c y 



unit (t ) for the 
a secondary 
estimated time 
to provide a 



(3) and for 

(4) and an 



o 1 after as few 



means (2) comprise 
tf e control unit (1 ) 
vith the following 
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- the update unit (5) is provided with means for calculating the update signal (upv(k)) as follows: 



up v (k) 



ijk)e(k) 
{a'*i„(k)r v (k) ) 



where: 

7 = step size 

a' - small constant to prevent division by 0; 
• the control unit (t) is provided with means for calculating the filter coefficients <9 CrV (fr)) as follows: 

System according to Claim 1, 2 or 3, characterised in that the system also comprises estimating means 
(6) for estimating the transfer function (H) of the secondary path, which means have an input connected 
to the sensor means (4) and an output connected to the adaptation unit (7) for transmitting an output 
signal for up-dating of the transfer function (H) used in the adaptation unit (7). 

System according to one of the preceding claims, characterised in that the control unit (1), the update 
unit (5) and the adaptation unit (7) are implemented together in the software of a computer. 

System according to Claim 4, characterised in that the control unit (1), the update means (5), the 
estimating means (6) and the adaptation unit (7) are implemented together in the software of a 
computer. 

System according to one of the preceding claims, characterised in that the cancellation-generating 
means (2) are selected from the following means: 

at least one loudspeaker, at least one vibration generator, at least one adjustable flow-induced 
noise source or a combination thereof, and in that the sensor means (4) are selected from the following 
means: 

at least one noise pressure sensor, at least one vibration recorder or a combination thereof. 
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